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FOREWORD 


This Final Report is submitted in fulfillment of the re- 
quirements of Article IV of NASA Contract NAS4-2812, Canadian 
Commercial Corporation Contract IPD. 70E5-80-1, SN 7PD80-00101. 
The work was conducted under the direction of the HiMAT Project 
Office, NASA/DFRC, Edwards AFB, California. The cognizant pro- 
pulsion engineer was Mrs. J.L, Baer-Riedhart , NASA/DFRC. 

Altitude facility data were provided from tests conducted at 
NASA Lewis Research Center. Mr. L.A. Burkardt was the NASA pro- 
ject engineer (analysis). 

The optimization of the thrust algorithm calibration was 
conducted by Computing Devices Company, Propulsion Systems 
Group. Mr. M.J. Hamer was the project engineer and performed 
the analysis in conjunction with Mr. R.I. Alexander. 
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INTRODUCTION 


BACKGROUND 

In-flight gross thrust computation is a continuing task for 
all agencies concerned with the development and test of air- 
craft. Computing techniques usually fall into the categories of 
direct force measurement, internal or gas generator calibration, 
and nozzle exit pressure traversing. The gross thrust computing 
technique developed by Computing Devices Company (ComDev) is a 
simplified internal one that uses a flow calibration of the en- 
gine tailpipe alone. 

SIMPLIFIED THRUST COMPUTING TECHNIQUE. The simplified tech- 
nique uses measured total and static pressures in the engine 
tailpipe and ambient static pressure to compute gross thrust. 
The equations are based on a one-dimensional analysis of the 
flow. The gas flow model accounts for friction, heat and mass 
transfer, and three-dimensional effects through the use of em- 
pirically-determined calibration coefficients. Instrumentation 
bias may also be eliminated by calibration. 

Gas generator methods for computing thrust require many en- 
gine measurements, comprehensive analytical or model work and 
extensive full-scale calibration testing. The simplified tech- 
nique requires fewer and simpler engine measurements and com- 
puting requirements are negligible in comparison to gas gener- 
ator methods. Thrust may be processed on-line in the aircraft 
without compromising accuracy. As a result, the simplified 
technique is considered to be suitable for flight test and for 
use on production engines. 

The technique was originally developed using ground level 
engine data on the J85-CAN-15 afterburning turbojet engine [1]. 
The method was later extended to two afterburning turbofan en- 
gines (TF30, FIDO) using NASA altitude facility data. The tur- 
bofan results [2j evaluated the technique over a wide range of 
Mach/altitude test conditions. The FlOO system accuracy was 
verified on a second engine without altering the calibration of 
the algorithm. The results were employed in-flight on a NASA 
F“15 test aircraft with an accuracy comparable to the engine 
manufacturer's gas generator method [3]. 
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HiMAT EVALUATION PROGRAM 


OBJECTIVE. The thrust computing system (TCS) for HiMAT is 
being jointly developed by NASA and ComDev with an objective to 
produce a system which will provide optimum accuracy for the 
computation of in-flight gross thrust of the HiMAT engine over 
the operational envelope during all steady-state modes of oper- 
ation. 

SCOPE. Under previous C4j NASA Contract NAS4-2644, ComDev 
provided engine pressure instrumentation design, manufacture and 
installation, system error analysis, and a preliminary gross 
thrust algorithm. Approximate values of the algorithm calibra- 
tion coefficients were obtained by using engine pressures and 
thrust values from the J85-GE-21 model specification, and pres- 
sure ratios from previous J85-CAN-15 experience. The prelimi- 
nary algorithm was used to estimate pressure transducer ranges 
and system sensitivity to pressure measurement errors prior to 
final calibration of the system using altitude facility testing. 

The current HiMAT work is discussed in this Report and uses 
altitude facility test data on one J85-GE-21 engine in order to 
optimize the gross thrust algorithm calibration coefficients. 
This work is considered significant for three major reasons: 1) 
it allows the afterburning turbojet algorithm to be evaluated 
over the engine envelope (previous J85 work was based on ground 
level data only), 2) the HiMAT engine was operated using differ- 
ent engine control schedules so that the simplified technique's 
accuracy can be evaluated with change in exhaust nozzle sched- 
ules (this work was recommended in [2]), and 3) the current 
HiMAT work investigated casing as well as liner static pressure 
taps. Previous ground level J85 tests used casing static pres- 
sure taps only. Finally, the current results are presented in 
such a way that the reader can assess the method's accuracy for 
a calibration based on data from one test condition versus the 
accuracy for a calibration based on several flight test condi- 
tions. The final gross thrust algorithm was delivered to NASA 
as a FORTRAN subroutine. 
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ENGINE TESTS 


ENGINE DESCRIPTION 

The HiMAT ground test engine is a J85-GE-21 afterburning 
turbojet, S/N 225326, equipped with a standard bill-of-materials 
(BOM) control system.. The engine's aft sections were modified 
for installation of the TCS hardware. The modifications were: 
1) installation of total pressure rakes in the center-cone sup- 
port body, 2) installation of static pressure taps in the after- 
burner casing and afterburner liner aft of the afterburner 
f lameholder s, 3) installation of static pressure probes at the 
nozzle inlet replacing the respective stand-off bolts, and 4) 
installation of harnesses and tubing for transferring the pneu- 
matic signals to an engine-mounted bulkhead at the compressor 
case. For engine operation requiring off-design exhaust nozzle 
scheduling, a separate throttle signal was sent to the after- 
burner fuel controller, which also controls the nozzle area, to 
command the nozzle to the desired area schedule. 


ALTITUDE TEST FACILITY 

A photograph of the J85-GE-21 engine installed in the NASA 
Lewis PSL3 altitude facility is shown in figure 1. The station 
locations and the instrumentation used in the facility are shown 
in figure 2. The facility had a calibrated load cell thrust bed 
for determining actual gross thrust. 

NASA provided estimates of bias error, precision error and 
total uncertainty for the facility measurement of gross thrust. 
The bias error is a constant 49N and was estimated from data 
system and instrumentation specifications. The precision error 
varies as a function of altitude and power setting and repre- 
sents the standard deviation of 40 consecutive data scans at a 
given test condition. Total uncertainty was calculated ac- 
cording to the method in [5]. For simplification the results 
were plotted as a function of thrust level as shown on figure 
3. The total uncertainty, through the precision error compo- 
nent, contains a contribution due to engine thrust fluctuations 
at a fixed power setting. As a result, this total uncertainty 
represents a conservative estimate of the facility thrust meas- 
urement accuracy. The actual facility accuracy lies between the 
bias limit and total uncertainty plotted on figure 3. 


ENGINE INSTRUMENTATION 

The engine instrumentation kit produced for the J85-GE-21 
HiMAT engine was adapted from a previous design used on J85-GE-5 
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and J85-CAN-15 engines which have accumulated in excess of 
10,000 flight hours. Total pressure is measured by four, 
three-probe rakes at the turbine exit (designated Pt 5 ) . Sta- 
tic pressure is measured in the afterburner flameholder region 
both on the liner (Ps 6L) and on the casing (Ps6C) each with 
4 pressure taps. Static pressure is measured on the liner at 
the nozzle entry region (Ps7) using 4 pressure taps which re- 
place liner stand-off bolts. Engine tailpipe measurement sta- 
tions are shown on figure 4. At each station the pressures are 
manifolded to provide a pneumatic average, and routed (using 
0.32 cm outside diameter tubing) to an outlet bracket mounted on 
the compressor case. The low profile of the pneumatic plumbing 
did not interfere with the installation of the engine in the 
vehicle . 

The detailed design of the pressure probes for Pt5, Ps6C 
and Ps7 remained unchanged from that used on the J85-GE-5 and 
J85-CAN-15 engines. For the HiMAT program, a new liner tap for 
station 6 was designed, ground tested, and used during the alti- 
tude testing. A photograph of the new liner tap design is shown 
in figure 5. The detailed design of all of the probes is given 
in [4]. 

The pressure probe and manifold axial and circumferential 
exact locations are shown on figure 6. Separate manifolds were 
used at station 6 to measure Ps6C and Ps6L* Ambient static 
pressure was determined from nozzle exit external static taps as 
shown. 

NASA provided estimates of bias error, precision error and 
total uncertainty for the pressure transducers used for meas- 
uring Pt5, PS6L/ Ps6C» Ps7 and PsO* The bias error is 
a constant 0.023 N/cm2 for Pt5» PS6L/ PS6C and Ps7 and 
a constant 0.007 N/cm2 for PsO* The bias errors were esti- 
mated from data system and instrumentation specifications. The 
precision error varies as a function of altitude and power set- 
ting and represents the standard deviation of 40 consecutive 
data scans at a given test condition. For simplification, the 
results were plotted as a function of pressure level as shown on 
figure 7. Since the total uncertainty contains a contribution 
due to engine pressure fluctuations at a fixed power setting, 
the actual pressure measurement accuracy lies between the bias 
limits and total uncertainty limits shown on figure 7. 

An additional set of pressure transducers for the TCS was 
supplied by NASA Dryden. The NASA Lewis transducers were found 
to be more repeatable and accurate during the altitude tests, 
therefore, they were used as the primary instrumentation. A 
comparison of the two sets of transducers was made and is shown 
in Appendix B. 
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TEST CONDITIONS 


The conditions for the standard nozzle test and open-sched- 
uled nozzle test are illustrated in figure 8. Table 1 lists 
these conditions including the power range for each test. All 
of the test conditions were at standard day temperature except 
for the Mach 0.4, 6100 m condition which was approximately 25® C 
hotter than standard. As shown, the engine was tested over a 
wide range of conditions including the extremes of the standard 
day engine envelope. The engine was tested with a clean inlet 
configuration (no distortion screens). The general test proce- 
dure was to establish the facility flow on a given Mach/altitude 
test condition. Once the engine was established (4 minutes) at 
each throttle setting, a data point was taken followed by a re- 
peat data point (l/2 minute later). 

The standard nozzle schedule and open nozzle schedule are 
shown on figure 9 for power lever angles from idle (0-13 degrees 
PLA) to military power (90-93 degrees PLA) . In the open mode, 
the engine exhaust gas temperature control at military power is 
downtrimmed by approximately 110°K and nozzle area is opened up 
by approximately 13%. At each test condition in the standard 
mode, typically 5 non-afterburning and 3 afterburning power set- 
tings were tested. In the open mode, only non-afterburning pow- 
er settings were tested. 
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DATA REDUCTION 


INTRODUCTION 

The engine batch data were received from NASA Lewis and were 
reduced to examine the consistency between the measured tailpipe 
ratios and between pressure ratio and facility gross thrust. 
This was done to identify any outliers in the data prior to 
calibration of the simplified gross thrust algorithm. 


DATA BASE 

Measured tailpipe pressures Pt5/ PS6L/ PS6C# PS7 and 
PSO and facility thrust Fqm were examined. This was done 
for all of the operating points# a total of 388 data points in- 
cluding the repeat points. 


BEHAVIOUR OF TAILPIPE PRESSURES 

Typical values of the J85-GE-21 tailpipe pressures are shown 
in figure 10 to show their general behaviour. The standard mode 
pressures are plotted as a function of power setting for the 
0.9, 7620 m test condition. The pressures increase in going 
from flight idle to a maximum near military power. Pt 5 and 
PS6L remain constant throughout afterburning. Ps6C drops 
slightly with degree of afterburning while Ps7 decreases mark- 
edly with degree of afterburning. 


TAILPIPE PRESSURE DIAGNOSTICS 

Liner pressure ratio Ps 6L/PT5 was plotted as a function 
of Ps 7/PT5 as shown on figure 11. The data were found to be 
very consistent, particularly in non-afterburning, over the full 
operating envelope. Only three data points out of the 388 were 
identified as measurement outliers, based on a three-sigma lim- 
it, and were eliminated. All three occurred at flight idle in 
the top left-hand portion of the envelope. As shown on figure 
11, the non-afterburning data points for both engine control 
modes collapse to a straight line with a spread of +0.57 per 
cent of the point at the 2-sigma level. The afterburning data 
also collapse to a straight line and the spread is +1.81 per 
cent of the point at the 2-sigma level. 

Casing static pressure ratio Ps6c/PT 5 was also plotted 
as a function of Ps 7/PT5 as shown on figure 12. These data 

were also found to be very consistent, with no measurement out- 
liers found. As shown on figure 12, the non-afterburning data 
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collapse to a straight line with a spread of +0.85 per cent of 
the point at the 2-sigma level (versus +0.57 per cent of the 
point for the liner data). The afterburning data, collapsed to 
within +1.27 per cent of the point (versus +1.81 per cent of the 
point for the liner data) . Compared to liner data, the casing 
data were less consistent in non-afterburning and more consis- 
tent in afterburning. 


THRUST VS PRESSURE DIAGNOSTIC 

To identify any outliers in the measured gross thrust data, 
facility corrected gross thrust, FQM/d# was plotted as a func- 
tion of pressure ratio Pts/PsO as shown on figure 13. No 
outliers in the measured thrust data were found. As shown on 
figure 13, the non-afterburning data points for both engine con- 
trol modes collapse to a straight line with a spread of +4.3 
per cent of the point at the 2-sigma level. Afterburning data 
could not be checked in this manner since a single correlation 
between thrust and Pts/PsO does not exist. 


ACCURACY PREDICTIONS 

Accuracy predictions in the next section were based on the 
bias error (b) and twice the precision error (2s) of the data. 
These errors were then combined using the method in [5] to pro- 
duce a total uncertainty (U) of the algorithm, which includes 
test stand uncertainty. 

U = |b| + 2s 
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ALGORITHM CALIBRATION RESULTS 


PINAL GROSS THRUST ALGORITHM 

The equations used in the final gross thrust algorithm are 
summarized below. Derivations of the basic equations are given 
in [1]. The algorithm is shown schematically in figure 14. 

ENGINE DATA REQUIRED. The following engine geometrical data 
are specified: 


CA7 = A 7 


CALIBRATION COEFFICIENTS. The following calibration coeffi- 
cients are determined by engine calibration, and their numerical 
values are supplied to the gross thrust algorithm: 

C56, A57, C 6 A, C 6 B, C 6 C, EA, EB 

An appropriate average value (y) for the ratio of specific heats 
of the exhaust gas is also defined during calibration. 

MEASUREMENTS. The following pressure measurements are re- 
quired within the engine tailpipe: 

Pt5» PS6L» PS7 

Ambient static pressure is also required: 

Pso 


If the engine is operating in a non-standard control mode, al- 
ternate values of A57, C 6 A, C 6 B, C 6 C, EA and EB may be selected 
using the KMODE signal. 

CALCULATIONS. The following calculations are required to 
determine the gross thrust. 

(a) Calibration coefficient E is calculated from: 

E = EA.Pgg + eb 


and is used to modify the measured value of Ps 7 to 
form a new value, Ps 7 ns 

Ps7n = f(PS7, E) 


8 


AOll/FR 



PT6 is calculated following the approach in [1]: 
PT6 = f(PT5/ PS6» C56, 7 ) 


The calibration procedure provided a relationship be- 
tween C67 and an estimated value of As/Ay, 
(Ag/A7)E. 

C67 = C 6 A.(Ag/A^)g + C 6 B.(Ag/A^)^ + C 6 C 


This was done to obviate the need for iteratively com- 
puting C67 as a function of the final computed value of 

^ 8 . (A 8 /A 7 )e is first determined from an esti- 
mated value of X 7 , X 7 E, using the isentropic flow 
relation : 


^^8/A7)e = f (X7E, y ) 


where X 7 E is determined using a combination of the 
isentropic flow relation: 


X 


7E 




y-1 

y 


1 


and the following equation, which uses an approximate 
value of C67 (= A57), to first calculate Pt 7 : 

Pt 7 = f(PT 6 » PS7n» ^57, 7 ) 


The approximations in this step are absorbed in the 
calibrated relationship for C67 as a function of 
<As/A7)e. 

Xe is calculated using the isentropic flow relation: 



X 7 is calculated using the exact value of C67 follow- 
ing the approach in [IJ: 
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X? = f(PS6, X6, PS7m C67) 


(f) Pt 7 is calculated using the isentropic flow relation: 


P 


T7 


^S7n 


1 + 


y-1 

2 y 



y 

y-1 


and Pt 8 = PT7 (assumed). This assumption is ab- 
sorbed by the calibration coefficients. 

(g) Pts/PsO is calculated to check for choking. For 

choked flow, defined by 



(y+i) 

2 


y 

y-1 


gross thrust is calculated using the final thrust equa- 
tion in ClD for choked flow: 

Fg = fl(Ps7n» X 7 , A 7 , PsO/ y ) 


(h) If the flow is unchoked, gross thrust is calculated 
using the final thrust equation in [ 1 ] for unchoked 
flow : 


Fg = f2(PS7n» X 7 , Ay, Pso» 7 ) 


The preliminary algorithm was calibrated using data at 
the Mach 0.9, 7620 m altitude test condition. This provided 
reasonable accuracy over the flight envelope. A significant im- 
provement in accuracy was made by recognizing that one of the 
coefficients (E) varied as a linear function of tailpipe pres- 
sure level. As a result, the preliminary algorithm was modified 
accordingly. Other revisions included the capability to use se- 
parate sets of calibration constants for each engine control 
mode (KMODE) and removing the requirement for having an after- 
burner check [ 6 ]. 


ALGORITHM CALIBRATION 

The HiMAT thrust algorithm (HIMATF) was calibrated on each 
of 7 sets of HiMAT J85-GE-21 engine pressure data and NASA Lewis 
measured thrust and ambient pressure. The seven sets of data 
were collected at the following conditions: 
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Mach = 0.9, 7620 m altitude, standard and open nozzle 
modes, liner Psg taps 

Mach = 0.9, 7620 m altitude, standard nozzle mode, 
liner Ps6 taps 


Mach = 0.9, 7620 m altitude, open nozzle mode, liner 
^S6 taps 


All Mach/altitudes, standard and open nozzle modes, 
liner Ps6 taps 


All Mach/altitudes, standard nozzle mode, liner 
Ps6 taps 

All Mach/altitudes, open nozzle mode, liner 
PS6 taps 


All Mach/altitudes, standard nozzle mode, casing 
Ps6 taps 


The calibration method produced a set of calibration coeffi- 
cients which enabled the thrust algorithm (HIMATF) to compute a 
value for thrust which was as close as possible to the measured 
thrust. The calibration coefficients were then placed in the 
thrust algorithm and thrust was computed using pressures from a 
set of data points different from those used to calibrate. This 
prediction of thrust was compared to the measured thrust. A 
comparison of the calibration coefficients E and C67 from the 
seven calibration data sets is shown in figure 15. Figure 16 
illustrates the scatter of the calibration coefficient C67 for 
the seven calibration sets; the scatter being proportional to 
the accuracy of the thrust algorithm. 


CALIBRATION AT ONE MACH/ ALTITUDE CONDITION 

The HiMAT thrust algorithm was calibrated at the Mach 0.9, 
7620 m altitude condition using both repeat and non-repeat data 
points. Points were deleted when the computed pressure ratio 
^TS/PSO was less than 1.65 for standard nozzle mode and less 
than 1.50 for open nozzle mode due to accuracy deterioration of 
the algorithm at low pressure ratios. Three calibrations were 
conducted at this Mach/altitude condition and the resulting 
calibration coefficients are listed in table 2. The first cali- 
bration (calibration #1) used both standard mode data at Mach 
0.9, 7620 m (runs lA, 5A, 9A) and open data (runs lOB and 14B) 
for a total of 59 points. This calibration was used to predict 
on all Mach/altitude conditions, both standard and open modes. 
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non-repeat points only. Duplicate runs lA, 5A and 14B at Mach 

0. 9, 7620 m and IB and 5B at Mach 0.6, 9140 m were deleted from 
the prediction to ensure equal representation of each Mach/alti- 
tude condition. Line 1 of table 4 shows the results of this 
prediction. With a total of 128 points, the average bias error 
of (FgC“Fgm)/PGM was -0.16 per cent and twice the pre- 
cision error was 2.27 per cent where Fqc is computed thrust, 
^GM is Lewis altitude facility measured thrust and precision 
error is the standard deviation of (Fqc-Fgm) /Fgm about the 
average bias. The bias and precision of Fqc-Fgm in kN are 
also given. The 128 predicted points are plotted in figure 17. 

Figure 18 and line 2 of table 4 show the results of predict- 
ing with calibration #1 on only the standard mode points in line 

1. The number of points decreased to 93, the bias error changed 
to -0.29 per cent and twice the precision error decreased to 
2.18 per cent. The prediction on open mode data only is plotted 
in figure 19 and is summarized in line 3 of table 4. For 35 
points, bias error was 0.21 per cent and twice the precision er- 
ror was 2.39 per cent. 

Calibration #2 was obtained by calibrating on only the 
standard mode points from Mach 0.9, 7620 m used in calibration 
#1. Line 4 of table 4 shows the prediction using calibration #2 
on the set of standard mode data only from all Mach/altitude 
conditions used in the line 2 prediction. The results of this 
prediction are shown in figure 20. Compared to line 2, bias er- 
ror remained about the same at -0.31 per cent but twice the pre- 
cision error improved to 1.83 per cent. When only open mode 
data from Mach 0.9, 7620 m were used to calibrate (calibration 
#3, line 5, table 4) there were only 17 data points available, 
possibly explaining the prediction results (figure 21) on the 35 
open mode data points from all Mach/altitudes. Compared to line 
3 of table 4, where both modes with a total of 59 points were 
used to calibrate, the bias error increased to 0.26 per cent and 
twice the precision error increased to 2.44 per cent. 


CALIBRATION AT ALL MACH/ALTITUDE CONDITIONS 

The HiMAT thrust algorithm was calibrated on data from all 
Mach/altitude conditions using only repeat points and deleting 
all points for which the computed pressure ratio Ptb/PsO was 
less than 1.65 for standard mode and less than 1.50 for open 
mode. The calibration coefficients for the three calibrations 
conducted at all Mach/altitudes are listed in table 2. The 
calibration coefficients for the calibration using the casing 
data are also shown and are discussed below. The first calibra- 
tion (calibration #4) used both standard and open mode data with 
runs lA, 5A, IB, 5B and 14B deleted. Points having Pt8/PS0 
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less than 1.65 for standard and less than 1.50 for open mode 
were also deleted leaving a total of 123 points. This calibra- 
tion was used to predict on a set of data similar in all re- 
spects to the calibration set with the exception that only 
non-repeat points were used. This prediction is described in 
line 1 of table 5 and the results are shown in figure 22. For a 
total of 128 predicted points# bias error was -0.16 per cent and 
twice the precision error was 2.10 per cent. Comparing this 
line to line 1 of table 4 shows that prediction accuracy im- 
proves slightly when more conditions are used to calibrate. 
Twice the precision error decreased from 2.27 per cent in table 
4 to 2.10 per cent in table 5. The bias error change between 
the two predictions was negligible. 

When calibration #4 was used to predict on only the standard 
mode data of line 1, table 5, the bias error changed to -0.20 
per cent and twice the precision error decreased to 1.97 per 
cent for a total of 93 points. Line 2 of table 5 and figure 23 
show the results of this prediction. The small reduction in 
precision error and bias error resulting from calibrating on 
several Mach/altitude conditions instead of only one is again 
evident when line 2 of table 5 is compared to line 2 of table 
4. Bias error decreased from -0.29 per cent to -0.20 per cent 
and twice the precision error decreased from 2.18 per cent to 
1.97 per cent. 

Predicting on only the open mode data in line 1 of table 5 
produced a bias error of -0.06 per cent and twice the precision 
error of 2.43 per cent for a total of 35 points. These results 
are illustrated in figure 24 and line 3 of table 5. 

The second calibration at all Mach/altitude conditions 
(calibration #5) used only standard mode data from line 1 of 
table 5. A total of 87 points were included in this calibra- 
tion. When this calibration was used to predict on the same 
standard mode data as in lines 1 and 2 of table 5, the resulting 
bias error was -0.20 per cent and twice the precision error was 
1.68 per cent. (See figure 25 and line 4 of table 5.) The im- 
provement in precision error resulting from specializing on en- 
gine operating mode is shown by comparing line 4 with line 2 of 
table 5. Twice the precision error decreased from 1.97 per cent 
to 1.68 per cent. Calibrating on several Mach/altitude condi- 
tions instead of one decreased twice the precision error slight- 
ly# on prediction, from 1.83 per cent to 1.68 per cent as shown 
in line 4 of table 4 and line 4 of table 5. 

The third calibration at all Mach/altitude conditions (cali- 
bration #6) used only open mode data from line 1 of table 5. 
Figure 26 and line 5 of table 5 show the results when calibra- 
tion #6 was used to predict on all the same open mode points as 
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in lines 1 and 3 of table 5. The bias error for the 35 points 
was 0.06 per cent and twice the precision error was 2.24 per 
cent . 

Tables 4 and 5 show that using data from several Mach/alti- 
tude conditions to calibrate results in slightly better predic- 
tions on all Mach/altitude data than if data from only one con- 
dition is used to calibrate. Calibrating on data from one con- 
dition/ both modes and predicting on all conditions, both modes 
produced a bias error of -0.16 per cent and twice the precision 
error of 2.27 per cent, whereas calibrating on all conditions, 
both modes and predicting on all conditions, both modes produced 
a bias error of -0.16 and twice the precision error of 2.10 per 
cent. Similarly, calibrating on only standard mode data at one 
condition and predicting on only standard mode data at all con- 
ditions produced a bias error of -0.31 per cent and twice the 
precision error of 1.83 per cent, while a calibration on all 
conditions produced a prediction bias error of -0.20 per cent 
and twice the precision error of 1.68 per cent. 

Optimum accuracy of the HiMAT thrust algorithm over the en- 
gine envelope is demonstrated in line 4 of table 5 for standard 
mode operation. The algorithm was calibrated using data from 
all Mach/altitude test conditions, standard mode only. Bias er- 
ror (b) was -0.20 per cent and twice the precision error (2s) 
was 1.68 per cent. These errors were combined using the method 
in [5] to produce an algorithm total uncertainty (U) of 1.88 per 
cent of the point, which includes altitude facility test stand 
uncertainty. 


U = |b| + 2s 


Line 5 of table 5 shows optimum algorithm accuracy for open 
mode operation. The algorithm was calibrated using data from 
all Mach/altitude test conditions, open mode only. Bias error 
was 0.06 per cent and twice the precision error was 2.24 per 
cent. Therefore, the total uncertainty of the algorithm in- 
cluding test stand uncertainty for open mode operation is 2.30 
per cent of the point. 


CALIBRATION ON OPEN MODE COMPARED TO NORMAL MODE 

Calibrations were produced for standard and open mode data 
together, only standard mode data and only open mode data. 
Calibrating on both modes using one test condition and predic- 
ting on standard mode data produced twice the precision error of 
2.18 per cent while calibrating on only standard mode data and 
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predicting on standard mode data produced twice the precision 
error of 1.83 per cent. Table 5 shows that when the calibra- 
tions used data from several Mach/altitude conditions the pre- 
dictions again were better when the operating mode was special- 
ized. Calibrating on only standard mode data and predicting on 
standard mode data produced twice the precision error of 1.68 
per cent instead of the 1.97 per cent resulting from the cali- 
bration on both modes. Similarly, an open mode calibration pre- 
dicting on open mode data produced twice the precision error of 
2.24 per cent instead of the 2.43 per cent resulting from the 
calibration on both modes. 


CASING TAP DATA COMPARED TO LINER TAP DATA 

At engine station 6, static pressure was measured both at 
the afterburner liner surface and in the space between the liner 
and the casing. Figure 27 compares casing with liner Ps6 at 
each test condition, standard mode. Figure 28 compares the same 
casing with liner Pgg data in a normalized manner. At maximum 
afterburning, casing Pg6 approximately equals liner Ps6 but 
at lower power settings, casing Ps6 is up to 7 per cent higher 
than liner Ps6. Figure 29 shows that in the open mode, casing 
^S6 varies up to 10 per cent above liner Ps6. 

Casing Ps6 pressures were used to calibrate the thrust al- 
gorithm at all Mach/altitude conditions, standard mode. The 
difference between liner and casing calibrations is evident in 
table 2 and in figure 15 in which the higher level and more 
negative slope of the casing C67 curve can be seen. The vari- 
able calibration coefficient E is also higher for casing than 
for liner pressures. The casing calibration was used to predict 
on casing pressure from all Mach/altitude, standard mode data. 
Figure 30 and line 6 of table 5 show that the prediction bias 
error was -0.26 per cent and twice the precision error was 3.10 
per cent for 93 points. Similar calibration and prediction 
(line 4, table 5) using liner Ps6 data produced twice the pre- 
cision error of only 1.68 per cent. 


THRUST ALGORITHM SENSITIVITY TO PRESSURE MEASUREMENT ERROR 

The sensitivity of the HiMAT thrust algorithm (HIMATF) to 
errors in the measurement of pressures Ps0» PS6» 

Pt5-Ps6, PT5~Ps7 (HiMAT flight transducer configuration) 
was calculated at several conditions. Sensitivities at the two 
power settings, military and maximum afterburning, were investi- 
gated for sea-level-static and six other Mach/altitude condi- 
tions covering the HiMAT engine operating envelope. 
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TRANSDUCER RANGES OF THE DRYDEN SYSTEM. Ambient static air 
pressure, Pso, is obtained from the aircraft data system and 
has a range of 0 to 13.8 N/cm2. Afterburner static pressure 
in the flameholder region, Ps6/ is measured by an absolute 
pressure transducer with a range of 0 to 34.5 N/cm2. The 
tailpipe pressure differential, Pt5-Ps 6» is measured by a 
differential pressure transducer with a range of 0 to 10.3 
N/cm2. The second differential pressure transducer, 

Pt5“Ps 7» also has a range of 0 to 10.3 N/cm2. Table 6 
lists the transducer ranges. 

THRUST ALGORITHM SENSITIVITY. Table 7 lists the 7 
Mach/altitude conditions and the change in computed thrust re- 
sulting from errors in the measurement of each pressure. The 
assumed pressure measurement errors are +0.25 per cent of full 
scale for Psq and +0.33 per cent of full scale for PS6/ 
Pt 5“PS6 and PT5“PS7‘ The thrust change was determined 
by computing thrust with the unperturbed pressures, and then 
perturbing one pressure by the amount of the measurement error 
and recomputing thrust. The change in thrust represents the 
thrust error resulting from the pressure measurement error for 
that particular pressure. This is done for all four pressures 
and the results are root-sum-squared to produce the estimated 
thrust change due to pressure measurement errors in all trans- 
ducers. 

Therefore at military power, standard day sea-level static 
conditions, pressure measurement errors of +0.25% of full scale 
for Pso and +0.33 per cent of full scale for Ps6/ 
Pt5~Ps 6 and PT5-PS7 produce a computed thrust change of 
+0.80 per cent of the point. At intermediate afterburning (PLA 
= 110®) the thrust change reduces to +0.71 per cent of the 
point. At Mach 0.9, 7620 m altitude, standard day, the military 

power thrust change is +1.06 per cent of the point and the maxi- 
mum afterburning change is +0.78 per cent of the point. The 
Mach/altitude condition at which the thrust algorithm is most 
sensitive to pressure measurement error is Mach 0.9 at 15240 m, 
where military power thrust change is +3.19 per cent of the 
point and maximum afterburning thrust change is +2.32 per cent 
of the point. 

The thrust changes resulting from pressure measurement er- 
rors are plotted in Figures 31, 32 and 33 for three Mach/alti- 

tude conditions for pressure measurement errors from 0 to 0.5 
per cent of full scale. 
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CONCLUSIONS 


1. The simplified gross thrust algorithm was evaluated for the 
J85-GE-21 engine using altitude facility data* Computed 
thrust values were compared to measured thrust values. The 
algorithm was found to be very accurate over the engine en- 
velope for both the standard engine mode and the open nozzle 
mode . 

2. The difference in the algorithm accuracy for a calibration 
based on data from one test condition is small compared to a 
calibration based on data from all of the test conditions. 

3. The algorithm accuracy is slightly improved when the cali- 
bration is optimized for each of the standard and open noz- 
zle engine operation modes. 

4. Greater accuracy was obtained using liner static tap data 
than using casing static tap data when optimum accuracy over 
the range of power settings is required. 

5. The algorithm based on the calibration set using all 
Mach/altitude test conditions for the standard mode produced 
the optimum accuracy (total uncertainty) over the engine en- 
velope for this mode which was 1.88 per cent of point. 

6. The optimum accuracy (total uncertainty) for the open nozzle 
mode was the result of the algorithm based on the all 
Mach/altitude open nozzle calibration set and was 2.30 per 
cent of point over the engine envelope. 
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RECOMMENDAT I ONS 


1. For cost effectiveness, the method could be calibrated using 
repeat runs from one Mach/altitude condition rather than the 
same quantity of data from a variety of Mach/altitude con- 
ditions . 

2. It was reported in [4] that the simplified approach has the 
potential to provide a +2 per cent of the point accuracy 
system over most of the flight test envelope. The altitude 
facility results confirm that with accurate pressure meas- 
urement (+0.33 per cent), the calibration and model error 
are small enough to produce the desired +2 per cent system. 
In pursuit of this goal, the following are areas for further 
investigation ; 

(a) Optimization of the algorithm accuracy on the HiMAT 
flight engines to account for any engine-to-engine 
biases. 

(b) Determination of the accuracy to include installataion 
effects. 

(c) Comparison of the algorithm based on ground level data 
only with the algorithm based on Mach/altitude data to 
predict thrust over the engine envelope. 

(d) Examination of the accuracy of the algorithm using 
casing tap data for ground level conditions on the 
three HiMAT engines. 

(e) Development of the HiMAT pressure measurement system to 
produce a throughput accuracy of +0.33 per cent, in- 
cluding transducers, data link and associated inter- 
faces. One method would be to control the transducer 
environment and use a double transmission scheme [4]. 
An alternate method would be to preprocess the pressure 
data before inputting to the data link, and to transmit 
signals proportional to thrust to the ground-based com- 
puter for final processing. These signals could be 
transmitted with less accuracy than the raw pressure 
data to achieve +2 per cent thrust accuracy. 
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Run No. 

Test condition 

Thrust 

Range 

(Throttle) 

Nozzle 

Mode 

Mach 

No. 

Altitude 

(m) 

lA' 

0.0 

sis 

Idle - 110 deg 

Standard 

lA 

0.9 

7620 

Idle - Max 



IB 

0.6 

9140 

Idle - Max 



2A 

0.9 

12190 

Idle - Max 



2B 

0.6 

6100 

Idle - Mil 



3A 

0.9 

9140 

Idle - Max 



3B 

0.9 

15240 

Idle - Max 



4A 

0.6 

10670 

Idle - Mil 



4B 

0.6 

3050 

Idle - Max 



5A 

0.9 

7620 

Idle - Max 



5B 

0.6 

9140 

Idle - Mil 



6A 

1.4 

12190 

90% - Max 



6B 

1.4 

13720 

90% - Max 



7A 

1.2 

10670 

90% - Max 



7B 

1.2 

13720 

90% - Max 



8A 

0.4 

3050 

Idle - Mil 



SB 

0.4 

9140 

Idle - 90% 



9A 

0.9 

7620 

Idle - Max 


f 

9B 

0.6 

9140 

Idle - Max 

Standard 

lOA 

0.9 

15240 

Idle - Mil 

Open 

lOB 

0.9 

7620 

Idle - Mil 



llA 

0.9 

9140 

Idle - Mil 



IIB 

0.6 

6100 

Idle - Mil 



12A 

0.6 

10670 

Idle - Mil 



12B 

0.4 

9140 

Idle - Mil 



13A 

0.6 

3050 

Idle - Mil 



13B 

0.4 

3050 

Idle - Mil 



14A 

0.9 

15240 

Idle - Mil 



14B 

0.9 

7620 

Idle - Mil 


f 

14C 

0.0 

sis 

Idle - Mil 

0, 

pen 

15 

0.4 

6100 

Mil 

Standard | 


Notes: !• 


lA* and 14C are “sea-level-static" tests or as close 
as PSL3 facility can achieve this condition. 


2. 15 is a reference condition used during tests. 


Table 1. Test conditions 
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Calibration No. 

1 

2 

3 

Mach/Altitude 

0.9/7620 m 

0.9/7620 m 

0.9/7620 ra 

Run No. 

lA, 5A, 9A, 
lOB, 14B 

lA, 5A, 9A 

lOB, 14B 

Liner or Casing 
PS6 

Liner 

Liner 

Liner 

Mode 

Standard 
and Open 

Standard 

Open 

Number of Points 

59 

42 

17 

GAMMA 

1.30 

1.30 

1.30 

CA7 

227.6 

1 

227.6 

227.6 

A57 

-0.19 

-0.16 

-0.29 

C56 

0. 570 

0.570 

0.570 

EA 

0.0011 

0.0011 

0.0011 

EB 

-0.0315 

-0.0165 

-0.0565 

C6A 

0.38289 

0.36943 

-0.42567 

C6B 

-0.48133 

-0. 58542 

0.63516 

C6C 

-0.07075 

0.06723 

-0.51846 


Table 2. Calibration coefficients (page 1 of 2) 
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ro 

to 


Calibration No* 

4 

5 

6 

7 

Mach/Alti tude 

All 

All 

All 

All 

Run No* 

All except 
lA, 5A, IB, 
5B, 14B 

All except lA, 
5A, IB, 5B 

All except 14B 

All except lA, 
5A, IB, 5B 

Liner or Casing 
PS6 

Liner 

Liner 

Liner 

Casing 

Mode 

Standard 
and Open 

Standard 

Open 

Standard 

Number of Points 

123 

87 

36 

87 

GAMMA 

1.30 

1*30 

1.30 

1.30 

CA7 

227.6 

227*6 

227.6 

227.6 

A5 7 

-0.19 

VO 

• 

O 

1 

-0.29 

0*20 

C56 

0.570 

0*570 

0.570 

0*570 

EA 

0.00085 

0.00085 

0.00110 

0*00085 

EB 

-0.02275 

-0.01275 

-0*05650 

0*02725 

C6A 

0.33337 

0.25997 

-0*66251 

2.41473 

C6B 

-0.46175 

-0.42871 

0*96826 

-5*03391 

C6C 

-0*04189 

0.01172 

-0*62773 

2*49395 


Table 2. Calibration coefficients (page 2 of 2) 




Run 

No. 

Mach 

No. 

Altitude 

(m) 

Nozzle 

Mode 

Plotting 

Symbol 

lA* 

0.0 

sis 

Standard 

A 

lA 

0.9 

7620 



O 

IB 

0.6 

9140 



tb. 

2A 

0.9 

12190 




2B 

0.6 

6100 



Q 

3A 

0.9 

9140 



O 

3B 

0.9 

15240 



A 

4A 

0.6 

10670 




4B 

0.6 

3050 



1 

5A 

0.9 

7620 



□ 

5B 

0.6 

9140 



o 

6A 

1.4 

12190 



0 

6B 

1.4 

13720 



o 

7A 

1.2 

10670 



m 

7B 

1.2 

13720 



V 

8A 

0.4 

3050 




8B 

0.4 

9140 

y 

r 


9A 

0.9 

7620 



O 

9B 

0.6 

9140 

Standard 

G 

lOA 

0.9 

15240 

Open 

A 

lOB 

0.9 

7620 



• 

llA 

0.9 

9140 



♦ 

IIB 

0.6 

6100 



m 

12A 

0.6 

10670 



♦ 

12B 

0.4 

9140 



m 

13A 

0.6 

3050 



I 

13B 

0.4 

3050 




14A 

0.9 

15240 



A 

14B 

0.9 

7620 

1 


■ 

14C 

0.0 

sis 

Open 

A 

15 

0.4 

6100 

Standard 

G 


Table 3. Plotting symbols 
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Accuracy 


Calibration Set 
(# of Data Points) 

Prediction Set 
(# of Data Points) 

Bias error, b 

Twice the 
precision 
error, 2s 




% 

(kN) 

% 

(kN) 


Calibration #1 
Mach 0* 9, 7620m. Standard 
and open modes. Use lA, 
5A, 9A, lOB, 14B. Repeat 
and non-repeat points. 

(59) 

All Mach/altitude test conditions. 
Standard and open modes. Non- 
repeat points* Use all data ex-* 
cept lA, 5A, IB, 5B, 14B. 

(128) 

-0.16 

(-0.006) 

2.27 

(0.195) 

2 

As in 1 

As in 1 but standard mode only* 

(93) 

-0.29 

(-0.016) 

2.18 

(0.210) 

3 

As in 1 

As in 1 but open mode only. 

(35) 

0.21 

(0.020) 

2.39 

(0.137) 

4 

Calibration #2 
As in 1 but standard 
mode only. 

(42) 

As in 1 but standard mode only. 

(93) 

-0.31 

(-0.012) 

1.83 

(0.179) 

5 

Calibration #3 

As in 1 but open mode only. 

(17) 

As in 1 but open mode only* 

(35) 

0.26 

(0.014) 

2.44 

(0.107) 


Table 4. Prediction accuracy using calibration data from one test condition 

(Mach 0.9, 7620n) 
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Accuracy 


Calibration Set 
(# of Data Points) 

Prediction Set 
(i of Data Points) 

Bias error# b 

Twice the 
precision 
error > 2s 




% 

(kN) 

% 

(kN) 


Calibration #4 
All Mach/altitude test 
conditions. Standard and 
open modes. Repeat points 
only. Use all data except 
lA, 5A, IB, 5B, 14B. 

(123) 

All Mach/altitude test conditions. 
Standard and open modes. Non-Re- 
peat points. Use all data except 
lA, 5A, IB, 5B, 14B. 

(128) 

-0.16 

(-0.014) 

2.10 

(0.172) 

2 

As in 1 

As in 1 but standard mode only. 

(93) 

-0.20 

(-0.017) 

1.97 

(0.185) 

3 

As in 1 

As in 1 but open mode only* 

(35) 

-0.06 

(-0.006) 

2.43 

(0.130) 

4 

Calibration #5 
As in 1 but standard 
mode only. 

(87) 

As in 1 but standard mode only. 

(93) 

-0.20 

(-0.014) 

1.68 

(0.158) 

5 

Calibration #6 

As in 1 but open mode only. 

(36) 

As in 1 but open mode only. 

(35) 

0.06 
( 0.000) 

2.24 

(0.095) 

6 

Calibration #7 
As in 1 but standard 
mode only. Casing Ps6* 

(87) 

As in 1 but standard mode only. 
Casing Ps6* 

(93) 

-0.26 

(-0.023) 

3.10 

(0.290) 


NJ 

cn 


Table 5. Prediction accuracy using calibration data from all Mach/altitude conditions 










Pressure 

Transducer 

Range 

(N/cm^) 

ComDev Specified 
Pressure Measurement 
Error (% F.S.) 

Pso 

0 - 13.8 

+0.25 

PS6 

0 - 34.5 

+0.33 

Pt5-PS6 

0 - 10.3 

+0.33 

PT5-PS7 

0 - 10.3 

+0.33 


Table 6. Pressure transducer ranges 


Power 

Setting 

Mach 

No. 

Altitude 

(m) 

Change 

in computed thrust (% of 

Point ) * 

Pso 

^S6 

p _p 

*^T5 ^S6 

^T5“**S7 

RSS 

Mil 

0 

0 

-0.182 

0.629 

-0.060 

0.456 

0.800 

110® PLA 

0 

0 

-0.188 

0.629 

0.021 

0.275 

0.712 

Mil 

0.9 

7620 

-0.216 

0.668 

-0.173 

0.781 

1.064 

Max AB 

0.9 

7620 

-0.237 

0.720 

0.075 

0.182 

0.783 

Mil 

0.6 

9140 

-0.358 

1.069 

-0.313 

1.282 

1.736 

Max AB 

0.6 

9140 

-0.383 

1.166 

0.128 

0.291 

1.268 

Mil 

0.9 

9140 

-0.256 

0.769 

-0.230 

0.949 

1.269 

Max AB 

0.9 

9140 

-0.281 

0.849 

0.086 

0.208 

0.922 

Mil 

0.9 

12190 

-0.386 

1.082 

-0.451 

1.654 

2.064 

Max AB 

0.9 

12190 

-0.419 

1.225 

0.086 

0.423 

1.365 

Mil 

1.4 

12190 

-0.227 

0.664 

-0.238 

0.907 

1.171 

Max AB 

1.4 

12190 

-0.249 

0.724 

0.060 

0.214 

0.797 

Mil 

0.9 

15240 

-0.654 

1.788 

-0.723 

2.455 

3.190 

Max AB 

0.9 

15240 

-0.718 

2.123 

0.152 

0.593 

2.323 


* Pressure measurement errors and ranges are shown in Table 6. 

Note : Sensitivity calculated for calibration #1, standard 

mode data. 


Table 7. Effect of pressure measurement errors on HiMAT 
thrust algorithm computed thrust 


26 


AOll/FR 




AOll/FR 





per cerrt 




Engine stations 


Engine pressufe piobes 

Pjg - Tuf bine discharge total pressure 

Flameholder static pressure 
Nozzle entry static pressure 


f’se - 

PcT - 


S7 


,S ,6 



Outlet bracket 
mounted on compressor 
front frame 


Figure 4 Tailpipe measurement stations 
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Figur« 7. NASA Lewis pressure tronsducer bios limit and total uncertainty 
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APPENDIX A 


PRESSURE PLOTS 


HiMAT J85-GE-21 engine tailpipe pressures, ambient static 
pressure and thrust were measured at several Mach/altitude con- 
ditions and power settings in both the standard nozzle mode and 
in the open nozzle mode at the NASA Lewis altitude test facil- 
ity. This Appendix presents plots of pressures and thrust as 
measured by the NASA Lewis altitude test facility transducers. 
The engine tailpipe pressures measured were Pt 5, PS6(Liner)» 
^S6( casing) and Ps7* Ambient pressure Pso and gross 
thrust Pqm were also measured. Data from all Mach/altitude 
test conditions and both standard and open modes are presented. 
Table A-1 lists the different pressure plots that are presented. 


Figure No. 

X 

Y 

A-1 

PS7/PT5 

PS6l/PT5 

A-2 

PS6L/PsO 

Fgm/5 

A- 3 

Pts/Pso 

Fgm/5 

A-4 

(PT5-PS6L)/PsO 

Fgm/^ 

A-5 

(PT5-Ps7)/Pso 

Fgm/ d 


Data are from NASA Lewis altitude facility tests on HiMAT 
J85-GE-21 engine using NASA Lewis pressure transducers, all 
Mach/altitude conditions, standard and open nozzle modes, 
all power settings. 


Table A-1. Pressure plots 
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APPENDIX B 


COMPARISON OF NASA DRYDEN TRANSDUCER RESULTS TO 
NASA LEWIS TRANSDUCER RESULTS 


During the HiMAT J85-GE-21 engine tests at the NASA Lewis 
altitude test facility, engine pressures Ps6» Pt5"Ps6 
Pt 5 -Ps 7 were measured by pressure transducers which were 
mounted near the engine and were designed for installation on 
the HiMAT vehicle. This NASA Dryden transducer package was 
independent of the NASA Lewis altitude facility pressure 
monitoring equipment. The laboratory-standard transducers of 
the Lewis altitude facility provided the pressure data which 
were used to calibrate the HiMAT thrust algorithm HIMATF. The 
Dryden pressure transducers also provided pressure data but were 
not used during the calibration since they were not 

environmentally controlled and their accuracy was therefore 
questionable. This Appendix compares Dryden and Lewis pressures. 

Figure B-1 compares liner Pgs measured by Dryden 

transducers to liner Pgg measured by the Lewis transducers for 
all Mach/altitude conditions, standard mode only, non-repeat 
points only and with runs lA, 5A, IB and 5B deleted. The 93 
points show an average bias of -0.15 per cent or -0.030 N/cra2 
and twice the precision error of 1.26 per cent or 0.147 N/cm2. 

Figure B-2 compares Dryden Pt5-Ps 6 to Lewis Pt5~PS6 

for the same points as in figure B-1. The average bias is 1.01 
per cent or 0.017 N/cm2 and twice the precision error is 2.91 
per cent or 6.964 N/cm2. Figure B-3 compares Dryden to Lewis 
Pt 5 -Ps 7 . The average bias is -1.21 per cent or -0.028 
N/cm2 and twice the precision error is 3.76 per cent or 0.087 
N/ cm2 . 

The HiMAT thrust algorithm, calibrated using Lewis pressures 
from the test condition at Mach 0.9, 7620 m altitude, standard 

and open mode, was used to predict thrust from Dryden pressures 
from all Mach/altitudes, both modes. Figure B-4 shows that 
prediction bias error was -0.88 per cent and twice the precision 
error was 3.66 per cent. This bias and increased precision 
error over the Lewis prediction is consistent with the bias and 
precision errors observed between Lewis and Dryden pressures. 
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